Dislocation core field. II. Screw dislocation in iron 
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FIG. 1. Screw dislocation periodic arrangements used for ab 
initio calculations: (a) T and (b) AT triangular arrangements, 
(c) quadrupolar arrangement, u and v are the unit cell pe- 
riodicity vectors, and A the dipole cut vector. In all cases, 



the Burgers vector b = ^ao [111] for easy, and 
hard core configuration. 
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The dislocation core field, which comes in addition to the Volterra elastic field, is studied for 
the (111) screw dislocation in alpha-iron. This core field, evidenced and characterized using ab 
initio calculations, corresponds to a biaxial dilatation, which we modeled within the anisotropic 
linear elasticity. We show that this core field needs to be considered when extracting quantitative 
information from atomistic simulations, such as dislocation core energies. Finally, we look at how 
dislocation properties are modified by this core field, by studying the interaction between two 
dislocations composing a dipole, as well as the interaction of a screw dislocation with a carbon 
atom. 

PACS numbers; 61.72.Lk, 61.72.Bb 



I. INTRODUCTION 

Ab initio calculations have revealed that a (111) screw 
dislocation in a-iron creates a core field in addition to 
the Volterra elastic field-. This core field corresponds to 
a pure dilatation in the {111} plane perpendicular to the 
dislocation line. It is responsible for a non negligible vol- 
ume change per unit length of dislocation line. The core 
field decays more rapidly than the Volterra field, as the 
displacement created by this core field varies as the in- 
verse of the distance to the dislocation line, whereas the 
displacement caused by the Volterra field varies as the 
logarithm of this distance. Such a dislocation core field 
is not specific to iron: a similar dilatation induced by the 
core of the screw dislocation can be deduced from the 
analysis of core structures obtained from first-principles 
in other body-centered cubic (bcc) metals, such as Mo 
and Ta^^''^. Atomistic simulations in other crystal struc- 
tures have also led to such a core field^"— . 

The purpose of this paper is to characterize this core 
field for a (111) screw dislocation in iron, and to see how 
it modifies the dislocation properties. In that purpose, 
ab initio calculations have been used to obtain the dislo- 
cation core structure. Then the dislocation core field has 
been modeled within the linear anisotropic elasticity the- 
ory, using the approach initially developed by Hirth and 
Lothe^, and generalized in the preceding paper^*^ to in- 
corporate the core field contribution into the dislocation 
elastic energy. This modeling allows extracting quantita- 
tive information from ab initio calculations, such as the 
dislocation core energy. Finally, the effect of the core field 
on the interaction of a screw dislocation with a carbon 
atom has been investigated, as well as on the properties 
of a screw dislocation dipole. 



II. ATOMISTIC SIMULATIONS 

A. Dislocation dipole 

Fully periodic boundary conditions have been selected 
to study the (111) screw dislocation in a-Fe. A dislo- 
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cation dipole is introduced into the simulation box, us- 
ing three different periodic distributions of dislocations. 
Within the triangular arrangement, initially proposed 
by Frederiksen and Jacobsen'^, the dislocations are po- 
sitioned on a honeycomb network [Figs. HJa) and [Ifb)] 
that strictly preserves the threefold symmetry of the bcc 
lattice along the [111] direction. Two different variants, 
which are linked by a 7r/3 rotation, are possible for this 
periodic arrangement: the twinning (T) [Fig. [Ha)] and 
the anti-twinning (AT) [Fig. [ijb)] triangular arrange- 
ment. The name of the variant refers to the fact that 
the dislocation dipole has been created by shearing a 
{112} plane either in the T or AT orientation. The third 
periodic arrangement, represented in Fig. [IJc), is equiv- 
alent to a rectangular array of quadrupoles. The peri- 
odicity vectors u and v defining these different arrange- 
ments are given in Ref.fll). For all periodic arrangements, 
the periodicity vector along the dislocation line is taken 
as the minimal allowed vector, i.e. the Burgers vector 
b = [111] where oq is the lattice parameter. Thus 
the number of atoms in the simulation box is directly 
proportional to the surface S of the unit cell perpendic- 
ular to the dislocation line. 

For each periodic arrangement, the dislocations are po- 
sitioned at the center of gravity of three neighboring [111] 
atomic columns. Depending on the sign of the Burgers 
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vector, two types of cores can be obtainedi^. In the easy 
core configuration, the hehcity of the lattice is locally 
reversed compared to the helicity of the perfect lattice. 
This configuration has been found to be the most stable 
in all atomistic simulations in bcc transition metals. In 
the hard core configuration, the three central neighbor- 
ing [111] atomic columns are shifted locally so that these 
atoms lie in the same (111) plane. This configuration 
actually corresponds to a local maximum of the energy 
of the dislocation arrangement. It can be nevertheless 
stabilized numerically by symmetry in the atomistic sim- 
ulations. The three periodic arrangements sketched in 
Fig- Hallow us to simulate a dipole whose both disloca- 
tions are either in their easy or hard core configurationii. 

The dislocation dipole is introduced into the simulation 
cell by applying the displacement field of each dislocation, 
as given by the anisotropic linear elasticity. A homoge- 
neous strain is also applied to the periodicity vectors, to 
minimize the elastic energy contained within the simula- 
tion boxiii^"— . This homogeneous strain corresponds to 
the plastic strain produced when the dislocation dipole 
is introduced into the simulation unit cell. It is given by 

_ fajAj + bjAj 
2S 

The orientations of the Burgers vector, b, and of the 
dipole cut vector. A, are defined in Fig. [1] Then 
the atomic positions are relaxed so as to minimize the 
energy of the simulation box computed by ab initio 
calculations^-. 

Two types of simulations have been performed in this 
work: simulations at constant volume and at zero stress. 
One can keep the periodicity vectors fixed and minimize 
the energy only with respect to the atomic positions. 
Within this constant volume simulation, the simulation 
box is subject to a homogeneous stress. Within the zero 
stress simulation, the unit cell is allowed to relax its size 
and shape, so that the homogeneous stress vanishes at 
the end of the relaxation. In both cases, the dislocation 
core field can be identified. 



B. Ab initio calculations 

The present ab initio calculations in bcc iron have 
been performed in the density functional theory (DFT) 
framework using the Siesta code^^, i.e. the pseudopo- 
tential approximation and localized basis sets, as in 
Refs. JJ, and iH. Comparison with plane-waves DFT 
calculationsiS has shown that this Siesta approach is 
reliable to study dislocations in bcc iron. The charge den- 
sity is represented on a real space grid with a grid spacing 
of 0.06 A that has been reduced after self-consistency to 
0.03 A. The Hermite-Gaussian smearing technique with 
a 0.3 eV width has been used for electronic density of 
state broadening. These calculations are spin-polarized 
and eight valence electrons are considered for iron. The 



Perdew-Burke-Ernzerhof (PBE) generalized gradient ap- 
proximation (GGA) scheme is used for exchange and cor- 
relation. A 3 X 3 X 16 fc-point grid is used for the dislo- 
cation calculations with unit cells containing up to 361 
atoms, and a 16 x 16 x 16 grid for the elastic constants. 

The obtained Fe lattice parameter is oq = 2.88 A, in 
good agreement with the experimental value (2.85 A). 
The DFT elastic constants are deduced from a fit on 
a fourth order polynomial over the energies for differ- 
ent strains ranging from —2 to 2%. This leads to the 
values of 248, 146 and 69GPa for the elastic constants 
Cu, Ci2 and C44 respectively, expressed in Voigt no- 
tation in the cubic axes. These values are close to the 
experimental ones, Cn — 243 and C12 — 145 GPa, ex- 
cept for the shear modulus C44 which is found stiffer 
experimentally (116 GPa). As a consequence, the elastic 
anisotropy within DFT is less pronounced than experi- 
mentally: DFT calculations lead to an anisotropic ratio 
A = 2C44/(Cii - C12) = 1.35 instead of 2.36. 

The three ab initio elastic constants yield to a shear 
modulus in {110} planes fiuo = (Cn — C12 + C44) /3 
equal to 57 GPa, which is close to the experimental value, 
Miio — 71 GPa. This parameter is of major importance 
for (111) dislocations as it controls their glide in {110} 
planes. Another important quantity is the logarithmic 
prefactor K = /x6^/47r controlling the main contribution 
to dislocation elastic energy. For a screw (111) disloca- 
tion, the shear modulus appearing in this prefactor is 
equal to 56 GPa for ab initio data and 64 GPa for experi- 
mental ones, thus in close agreement. The error between 
ab initio and experimental elastic constants should there- 
fore not affect too much our results. For consistency, all 
elastic calculations below are performed using ab initio 
elastic constants. 



III. CORE FIELD CHARACTERIZATION 

The dislocation core field can be modeled by an equi- 
librium distribution of line forces2S parallel to the dis- 
location and located close to its core^. For the (111) 
screw dislocation in iron, the center of this distribution 
corresponds exactly to the position of the dislocation, i.e. 
to the center of gravity of three (111) neighboring atomic 
columns, for symmetry reasons. At long range, and at a 
point defined by its cylindrical coordinates, r and 9, this 
distribution generates an elastic displacement given by a 
Laurent series (see preceding paper—) 

00 ^ 

n=l 

The main contribution of this series, i.e. the term n — 1, 
is completely controlled by the first moments Mij of the 
line force distribution. Knowing this second rank tensor 
Mij, one can not only predict the elastic displacement 
and stress associated with the core field^, but also the 
contribution of the core field to the elastic energy and to 
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FIG. 2. Homogeneous stress observed in ab initio calculations 
for the three periodic dislocation arrangements in the easy 
core configuration: (a) T and AT triangular; (b) quadrupolar 
dislocation arrangements. 

the dislocation interaction with an external stress field—. 
It is thus important to know the value of the first moment 
tensor Mij , and we will see how it can be deduced from 
ab initio calculations. 



A. Simulations with fixed periodicity vectors 

A homogeneous stress is observed in the ab initio cal- 
culations, when the periodicity vectors are kept fixed, 
and when only the atomic positions are relaxed. The six 
components of the corresponding stress tensor are shown 
in Fig. [5] for the three periodic dislocation arrangements, 
in the easy core configuration. The stress components 
are expressed in the axes = [112], By — [110] and 
— [111]. The main components of the stress tensor 
are axx and ayy, and the other components can be ne- 
glected. The stress components axx and ayy vary roughly 
linearly with the inverse of the number of atoms, and so 
with the inverse of the surface S of the simulation box. 
For the two variants of the triangular arrangement, T and 
AT, the stress components axx and ayy are exactly equal 



FIG. 3. Moments Mxx and Myy of the line force dipoles for 
the three periodic dislocation arrangements in the (a) easy 
and (b) hard core configuration. 

[Fig. [2la)]. Indeed, the threefold symmetry along the 
[111] direction obeyed by this dislocation arrangement is 
also imposed to the homogeneous stress. On the other 
hand, the quadrupolar arrangement breaks this symme- 
try. As a consequence, axx and ayy slightly differ from 
each other in this case [Fig. [5l^b)]. 

The two core fields of the dislocations composing the 
dipole are responsible for this homogeneous stress, as 
shown in Ref. U. If My is the first moment tensor of 
the line force distribution representative of the core field, 
the homogeneous stress in the simulation box is given by^ 

The factor of 2 in this equation arises from the fact that 
two dislocations constituting the dipole are introduced 
within the simulation box. 

As the stress components other than axx and ayy 
can be neglected, the force moment tensor M will have 
only two non-zero components, Mxx and Myy, as shown 
in Eq. ([T]). These two components deduced from the 
stress computed from DFT calculations are represented 
in Fig. Elja) for the easy core configuration, and in 
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Fig. [3jb) for the hard core configuration. Within the 
two triangular arrangements, the core field is a pure bi- 
axial dilatation {M^x = Myy), whereas the core field has 
a small distortion component (M^x — Myy ^ 0) within 
the quadrupolar arrangement. This distortion is asso- 
ciated with the broken threefold symmetry within the 
quadrupolar periodic arrangement. It may arise from a 
polarizability^^"^'^ of the core field: the moments AI char- 
acterizing the core field may depend on the stress applied 
to the dislocation core. Such a polarizability may also be 
the reason why the moments obtained within the T vari- 
ant of the triangular arrangement slightly differ from the 
moments obtained within the AT variant (Fig. [S]). 

One can also observe a dependence of the moments 
with the size of the simulation box. The cell size that 
can be reached in DFT does not allow us to obtain well- 
converged values for the moments. Nevertheless, the core 
field of the screw dislocation in iron can be considered as 
a pure biaxial dilatation of amplitude Mxx = Myy = 
650 ±50GPaA2. Interestingly, the easy and the hard 
core configurations are characterized by the same mo- 
ments [Figs. [SJa) and[3Jb)], although the atomic struc- 
tures of the dislocation cores are completely different be- 
tween these two configurations. This probably indicates 
that the dislocation core field does not arise from pertur- 
bations due to the atomic nature of the core, but rather 
from the anharmonicity of the elastic behavior. 



six components of the homogeneous strain are given by 
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B. Simulations with relaxed periodicity vectors 



where the elastic constants C, _ _ 
notation, in the axes Gx = [112], = [110], and 
Gz — [111]- This shows that, for positive moments Mxx 
and Myy, a dilatation perpendicular to the dislocation 
and a contraction parallel to the dislocation line are pro- 
duced. This exactly corresponds to what we observe in 
the ah initio calculations. Thus one can define a dislo- 
cation formation volume perpendicular to the dislocation 
line, 5V±_ — (e°2.-|-e°y)S'/2 and a formation volume paral- 
lel to the dislocation line, 5V\\ = e°^S'/2, where the values 
of the formation volume are defined per unit of disloca- 
tion line. The DFT results are given in Table ID The 
expressions of the homogeneous strain [Eq. ([2])] also lead 
to a ratio between the two dislocation formation volumes 
depending only on the elastic constants as 



Ab initio simulations, in which the periodicity vectors 
are allowed to relax so as to minimize the energy, have 
also been performed. In this homogeneous strain, 

e^j, can be computed. This strain is related to the core 
field of the two dislocations within the simulation box, 
through the relation^ 
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where the elastic compliances Sijki are the inverse of the 
elastic constants^'*. If the dislocation core field is as- 
sumed to be an elliptical line source of expansion charac- 
terized by the two non-zero moments Mxx and Myy, the 
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Using the values of the elastic constants calculated in 
DFT, we predict 5Vj_/(5V|| = —2.0. This is in reason- 
ably good agreement with the values obtained from the 
homogeneous strain computed in ab initio calculations 
(Table U. 

The moments Mxx and Myy characterizing the dislo- 
cation core field can be deduced from the homogeneous 
strain computed in DFT, using the system of equations 
([2). We choose to derive Mxx and Myy from the compo- 
nents Exx and Syy of the strain, and the resulting values 
are given in Table U for the three arrangements. These 
values are in good agreement with those derived from 
atomistic simulations with fixed periodicity vectors (Sec. 
IIII A[) . For all periodic arrangements of dislocations, the 
dislocation core field can be considered as a pure biax- 
ial dilatation, and we neglect the difference between Mxx 
and Af,,,,. 
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TABLE I. Formation volumes SV± perpendicular and SV\i par- 
allel to the dislocation line per unit of dislocation line com- 
puted in DFT for the three dislocation arrangements. A'^ is the 
number of atoms contained within the simulation box. M^^ 
and Myy are the moments of the dislocation core field deduced 
from the components ej^, and Syy of the homogeneous strain. 
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(b)Hard core configuration 
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IV. CORE FIELD IN ATOMISTIC 
SIMULATIONS 

Now that the dislocation core field has been character- 
ized, we examine its infiuence on atomistic simulations. 
First we look at the atomic displacements observed in ab 
initio calculations: part of this displacement arises from 
the core field. Then we show that the dislocation core en- 
ergies can be extracted from these calculations when the 
core field contribution is considered in the elastic energy. 



A. Atomic displacement 

Relaxation of atomic positions in ab initio calculations 
leads to the definition of atomic displacements induced 
by the periodic array of dislocation dipole. The displace- 
ment along the [111] direction, i.e. along the dislocation 
line, also called the screw component, is very close to the 
anisotropic elastic solution corresponding to the Volterra 
field. There is no substantial contribution of the core 
field on this displacement component. 

A displacement perpendicular to the dislocation line 
can also be evidenced in atomistic simulations, i.e. an 
edge component (Fig.H]). Part of this displacement com- 
ponent corresponds to the dislocation Volterra field and 
arises from elastic anisotropy. The dislocation core field 
also contributes to the edge component. The Volterra 
contribution is more long-ranged than the core field con- 
tribution, as the former varies with the logarithm of the 
distance to the dislocation, whereas the latter varies with 
the inverse of this distance. Nevertheless, both contri- 
butions evidence a similar amplitude in the simulations. 



because of the reduced cell size used in DFT calculations. 

We can subtract from the atomic displacement given 
by ab initio calculations, the displacement correspond- 
ing to the superposition of the Volterra and the dislo- 
cation core fields, as predicted by the anisotropic linear 
elasticity taking full account of the periodic boundary 
conditions'^. The resulting map, represented in Fig. |4j 
shows that the elastic modeling manages to reproduce 
the ab initio displacements for all periodic arrangements 
even close to the dislocation core. The anisotropic linear 
elasticity fails to reproduce the ab initio atomic displace- 
ment only on the atoms within the core: displacements in 
the core are too large for applying a small perturbation 
theory such as elasticity. 

B. Dislocation core energy 

The screw dislocation core energy is deduced from 
DFT simulations by subtracting the elastic energy to the 
excess energy of the unit cell computed by ab initio tech- 
niques. The elastic energy is calculated by taking into 
account the elastic anisotropy and the periodic bound- 
ary conditionsiS. 

If only the Volterra contribution is considered in the 
elastic energy, the resulting core energies strongly depend 
on the periodic arrangement (Fig. [5]). As shown in Ref. [H, 
one can not conclude on the relative stability of the two 
different core configurations of the screw dislocation. The 
AT triangular geometry and the quadrupolar geometry 
predict that the easy core is the most stable configura- 
tion, whereas the T triangular arrangement leads to the 
opposite conclusion. Then if both the Volterra and the 
core fields are considered in the elastic energy, the core 
energies do not depend anymore on the periodic arrange- 
ment. A cell size dependence has been evidenced (Fig. [5]), 
and in all geometries the easy core is more stable than the 
hard core configuration. The convergence is reached for 
a reasonable number of atoms: E^°'^'^ = 219 ± 1 meV A^' 
for the easy core configuration, and 227 ± 1 meV A"' for 
the hard core configuration. These core energies are given 
for a core radius Tc = 3 A, as this value of Tc has been 
found to lead to a reasonable convergence of the core 
energy with respect to the cell size (c/. Appendix [A)) . 

To better understand how the core energy converges, 
one can decompose the elastic energy into the differ- 
ent contributions, and look how these contributions vary 
with the length scale. If one neglects the dislocation core 
field, the elastic energy of the simulation box containing 
a dislocation dipole is given by 

^clas ^ 2Sclas ^ ^^^0 ^ ^into ^ (3) 

where K'^ is a second rank tensor, which depends only on 
the elastic constants, and Vc is the dislocation core cutoff. 
The two first terms on the right hand side define the 
elastic energy of the dipole contained in the simulation 
box: E°^^ is the core traction contribution-^^, and the 
second term corresponds to the cut contribution. iJy-'v 
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(a)T triangular periodic arrangement (361 atoms). 
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(b)AT triangular periodic arrangement (289 atoms). 
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(c)Quadrupolar periodic arrangement (209 atoms). 



FIG. 4. Planar displacement map of a periodic unit cell containing a screw dislocation dipolo in the easy core configuration, 
as obtained from ah initio calculations (left), and after subtraction of the Volterra and the core elastic fields (right). Vectors 
correspond to the (111) in-plane displacement and have been magnified by a factor of 50. Displacements smaller than 0.01 A 
are omitted. For clarity, displacements of the six atoms belonging to the dislocation cores are not shown on right panel. 
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(a) Easy core configuration 
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FIG. 5. Core energy of the screw dislocation in the (a) easy 
and the (b) hard core configuration. Sohd symbols correspond 
to the core energies obtained when only the Volterra field 
is considered, and open symbols, to the core energies when 
both the Volterra and the core fields are taken into account 
(rc = 3A). 



corresponds to the interaction of the dislocation dipole 
with its periodic images^^. If the core field is taken into 
accounijiS, the elastic energy becomes 



iclas 



clas 



Tc 



o cnintcr 



K-J. (4) 



where the fourth rank tensor , which only depends on 
the elastic constants, enables one to calculate the core 
field contribution to the dislocation self energy. iJy^" 
and E™}" correspond to the interaction of the disloca- 
tion core field with respectively the Volterra field, and the 
core field of the other dislocations, i.e. the second dislo- 
cation composing the dipole, as well as the image dislo- 
cations due to the periodic boundary conditions. When 
the periodicity vectors, u and v, and the dipole cut. A, 
are scaled by the same factor A, -E™-" varies as 1/A, and 
E'c'^c'^, as The number N of atoms in the sim- 

ulation box is proportional to A~. The comparison of 
equations Q and Q shows then that the neglect of the 



core field leads to a core energy that converges as N~^/^. 
In the case of dislocation periodic arrangements, which 
are centrosymmetric like the quadrupolar arrangement of 
Fig- W[c) , i?v-™ vanishes in Eq. ^ . The core energy con- 
verges thus as when the core field is neglected. In 
all cases, this convergence is too slow to extract a mean- 
ingful core energy from ab initio calculations. Moreover, 
the obtained value does not really correspond to the core 
energy E'^°'^° : it results from the summation of E'^°^''^ and 
the core field self elastic energy, MijKf-^.iMkil/rc^ . 

Another interesting feature comes from the linear de- 
pendency of Ey*^^ with the Burgers vector. Going from 
the easy to the hard core configuration of the dislocation 
dipole, i.e. inverting the sign of the Burgers vector, one 
only changes the sign of Ey*f^, whereas all other con- 
tributions to the elastic energy remain constant. The 
easy and the hard core configurations have indeed a core 
field with the same amplitude (Fig. [3]). The contribu- 
tion -By-^c positive for a dislocation dipole in its easy 
core configuration within the T triangular arrangement 
[Fig.[TJa)]. Therefore, when the core field contribution is 
not included into the elastic energy, one underestimates 
the stability of the easy core configuration with respect 
to the hard core configuration within this geometry. The 
AT triangular arrangement leads to the opposite conclu- 
sion [Fig.[ljb)], since becomes negative for the easy 
core configuration. 

The underestimation or overestimation of the stability 
of the easy core illustrates the importance of consider- 
ing the dislocation core field in the elastic energy, when 
extracting quantitative properties from atomistic simu- 
lations. This is especially true for ab initio calculations, 
in which the small cell size makes it difficult to obtain 
converged values. Such a conclusion is not restricted to 
the calculation of dislocation core energies. For instance, 
extraction from atomistic simulations of the Peierls en- 
ergy barriers, and of the associated Peierls stresses, will 
also require the complete modeling of the dislocation core 
field. 



C. Dislocation line energy and line tension 

It is interesting to evaluate the different contributions 
to the dislocation line energy. Using the developed elas- 
tic model, the energy^ of a straight screw dislocation 
contained in a cylinder of radius i?oo is 



E' 



Ticlas 



I 



-M,,Kf^,iMM — . (5) 

The core energy has been found to be E'^°''° = 219 ± 
ImeVA"""^ for the easy core configuration, and 227 ± 
1 meV A^"'^ for the hard core configuration, with Tc = 3 A. 
The Volterra elastic field leads to two energy contribu- 
tions: the contribution of the core traction^^ is Ef'^^ = 
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FIG. 6. Elastic contributions [Eq. ((5])] to the dislocation 
(a) line energy E{e) and (b) line tension T{e) = E{e) + 
E(6) /A9^ as a function of the dislocation character 6 in 
the {110} glide plane. 6 = Q for the screw orientation. 



1 meV A~^, for a dislocation cut corresponding to a {110} 
glide plane, and the cut contribution, which corresponds 
to the third term in Eq. (O, is equal to 1.60eVA~^, 
where we have assumed the ratio Rqo/tc = 10*, which 
corresponds to a characteristic dislocation density of 
10^^ m~^. Finally, the core field contribution to the elas- 
tic energy, which corresponds to the last term in Eq. ([5]) , 
is equal to 42meVA~^. It is clear that most of the dis- 
location energy arises from the Volterra elastic field, and 
is associated with the cut contribution. Other contribu- 
tions, which are all associated with the dislocation core, 
account for about 14 % of the dislocation energy. In par- 
ticular, the contribution of the core field is less than 3 %. 

The dislocation line tension is actually more important 
than the line energy as it controls the shape of dislocation 
loops and curved dislocationa2Z. It is defined as 



T{e) = E{9) 



d6'2 ■ 



if E{d) is the dislocation line energy [Eq. ([5])] written 
as a function of the dislocation character 9, i.e. the dis- 
location orientation. We can evaluate this line tension 
by considering all elastic contributions entering in the 
dislocation line energy and neglecting the dependence 
of the dislocation core energy with 9 which we do not 
know. For the dislocation core field, we assume that the 
dipole tensor My does not depend on the dislocation 
orientation as we do not have any information on such 
a variation: the variation with the dislocation orienta- 
tion of the associated line energy only arises from elastic 
anisotropy. The different elastic contributions to the dis- 
location line energy E{9) and the associated line tension 
T{9) are shown in Fig.[6]for a dislocation character going 
from an edge orientation {9 — ±90°) to a screw orienta- 
tion (9 = 0). The most important contribution to the line 
tension arises, once again, from the cut contribution of 



the Volterra elastic field. In view of the values obtained, 
it looks reasonable to neglect other elastic contributions, 
as usually done in line tension models^. Nevertheless, 
some other studies have obtained higher relative contri- 
butions of the core field'*'^, which may therefore have a 
more important effect on the line tension. 



DISLOCATION INTERACTION WITH A 
CARBON ATOM 



We examine in this part how the dislocation core 
field influences the dislocation properties by modifying 
the way a dislocation can interact with its environment. 
First, we study the interaction between a (111) screw 
dislocation and a carbon atom in a-iron. 

In Ref. [H, the binding energy between a carbon atom 
and a screw dislocation in iron as predicted by the linear 
anisotropic elasticity has been compared to the one cal- 
culated by atomistic simulations based on an empirical 
potential approaches,. A quantitative agreement between 
both modeling techniques was obtained as long as the C 
atom was located at a distance greater than 2 A from the 
dislocation core. This result evidences the ability of lin- 
ear elasticity to predict this interaction. The empirical 
potential for iron^i^ used in Ref. [H does not lead to 
any core field for the screw dislocation, at variance with 
the present ab initio results. It is interesting to include 
now the core field contribution into the elastic field of 
the screw dislocation, to investigate its influence on the 
interaction between a carbon atom and a screw disloca- 
tion. Such a contribution to the interaction energy be- 
tween a solute atom and a dislocation has already been 
considered in the case of a substitutional impurity by 
Fleischer— who showed that it partly contributes to the 
solid solution hardening. 



A. Carbon atom description 

First we need to deduce from ab initio calculations a 
quantitative representation of a carbon atom embedded 
in an iron matrix. A solute atom is modeled in elasticity 
theory by its dipolar tensor, Pij , which corresponds to the 
flrst moments of the equilibrated point force distribution 
equivalent to the impurity. This tensor is deduced from 
atomistic simulations of one solute atom embedded in the 
solvent (c/. Appendix iBjl . 

Carbon atoms are found in the octahedral interstitial 
sites of the bcc lattice. The dipolar tensor, Pij, expressed 
in the cubic axes, is diagonal with only two independent 
components, because of the tetragonal symmetry of the 
octahedral site. Three variants can be obtained, depend- 
ing on the orientation of the tetragonal symmetry axis. 
For the [001] variant of the C atom, ab initio calculations 
lead to Pxx — Pyy — 8.9 and P^z = 17.5 eV (c/. Appendix 



11). 



9 



^ 0.4 




FIG. 7. Binding energy between a 1/2[111](110) screw dis- 
location and a carbon atom for different positions x of the 
dislocation in its (110) glide plane. The C atom lies in a [100] 
octahedral site in the plane h — duo ~ 2.04 A above the glide 
plane. The binding energy is calculated using the anisotropic 
elasticity theory considering only the Volterra field, or both 
the Volterra and the core fields of the dislocation. 



B. Binding energy 

Linear elasticity theory predicts that the binding en- 
ergy between the carbon atom characterized by its dipo- 
lar tensor, Py, and the screw dislocation is given by 



iid p c"^ 



(6) 



where efj is the elastic strain created by the dislocation. 
We use the linear anisotropic elasticity to calculate sfj 
by taking into account only the Volterra field, or both 
the Volterra and the core fields created by the disloca- 
tion. We consider that the core field is created by the 
line force moments M^^ = Myy = 650 GPa previously 
deduced. In Fig. [71 we represent the variation of the 
binding energy when the dislocation glides in a {110} 
plane, while the carbon atom remains at a fixed distance 
h from the glide plane. The first derivative of the plotted 
function gives the force exerted by the C atom on the 
gliding dislocation. When the C atom is close enough 
to the dislocation, the core field modifies the binding en- 
ergy. In particular, the binding of the C atom is stronger 
in the attractive region when the dislocation core field is 
considered. Thus the pinning of the screw dislocation by 
the C atom is enhanced by its core field. Conversely, the 
dislocation core field leads to a stronger repulsion in the 
repulsive region. 

When the separation distance between the C atom and 
the screw dislocation is high enough (> 20 A), the dis- 
location core field does not affect anymore the binding 
energy. One can consider that the C atom interacts only 
with the dislocation Volterra field. 



[110] 




[111] [112] 



FIG. 8. Screw dislocation dipole. F is the force exerted by 
one dislocation on the other. 




FIG. 9. Variation of the elastic energy, Ai5, and of the x 
component of the force F acting on the dislocations, Fx, with 
the distance between both dislocations, x, for a screw dipole 
of height /i = 10 A. 



VI. PASSING PROPERTIES OF A SCREW 
DISLOCATION DIPOLE 

We look in this part at how the dislocation core field 
modifies the equilibrium properties of a screw dislocation 
dipole. Dislocation dipoles play a significant role in single 
slip straining, where they can control the material flow 
stress. Such a situation arises for instance in fatigued 
metals, where dislocations are constrained to glide in the 
channels between dislocation walls^^iM. The saturation 
stress of the persistent slip bands is then partly controlled 
by the critical stress needed to destroy the dislocation 
dipoles. 

We consider a screw dislocation dipole in bcc iron. The 
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dipole is characterized by the height between each dislo- 
cation ghde plane, h, and by the projection of the dipole 
vector on the glide plane, sketched in Fig. [H Then 
we calculate the variation of the interaction energy be- 
tween the two dislocations composing the dipole, A.E, 
when dislocations glide, i.e. h is kept fixed while x varies. 
This dislocation interaction energy is computed using 
the linear anisotropic elasticity's, and considering that 
the two dislocations composing the dipole interact only 
through the Volterra field, or through both the Volterra 
and the core fields. This variation of energy, AE, is repre- 
sented in Fig. [5] for a dipole height h= 10 A. The dipole 
equilibrium angle corresponds to the minimum of A.E, 
i.e. 9 = 0. This value predicted by the anisotropic elas- 
ticity is equal to the one given by the isotropic elasticity. 
This contrasts with what is found in fee metals, where 
the dipole equilibrium angle strongly depends on elastic 
anisotropy for a screw dislocation dipole"^^. The dislo- 
cation core field does not modify the dipole equilibrium 
angle (Fig. [5]). Nevertheless, when both the Volterra and 
the core fields are included into the computation of the 
dipole elastic energy, the energy that defines the dipole 
equilibrium becomes steeper than when the dislocation 
core field is omitted. Thus the attraction between both 
dislocations is stronger when the core field is taken into 
account. This is obvious when looking at the glide com- 
ponent of the force exerted by one dislocation on the 
other one, (Fig. [9]). This force is the first derivative of 
AE with respect to x. When the dislocation core field is 
included in the dislocation interaction, this force goes to a 
higher maximum value than when only the Volterra elas- 
tic field is considered. To destroy the dislocation dipole, 
one needs to exert on a dislocation a force that is higher 
than the force arising from the interaction with the other 
dislocation. This shows that the dislocation core field 
leads to a more stable dipole. 

If a homogeneous stress is applied, the gliding force on 
each dislocation is simply the Peach- Koehler force, hayz- 
As the applied stress is homogeneous, no force originates 
from the dislocation core fieldiS. Therefore the dipole 
passing stress, i.e. the applied stress needed to destroy 
the dislocation dipole, is given by the maximum of the 
glide component of the Peach-Koehler force divided 
by the norm of the Burgers vector. This passing stress 
depends on the dipole height, h, as shown in Fig.tTOl The 
inclusion of the core field into the dislocation interaction 
leads to a higher passing stress, especially for small dipole 
heights. But the eS'ect is relevant at a spacing where the 
dipole would certainly have cross-slipped to annihilation. 
For large dipole heights {h > 20 A), the dislocation core 
field does not influence too much the passing stress, and 
one can consider that dislocations interact only through 
the Volterra elastic field to calculate the passing stress. 

Without the dislocation core field, the dipole pass- 
ing angle, 6'inax, does not depend on the dipole height 
(Fig. [T0|) . The value given by the anisotropic elasticity is 
close to the 7r/4 value predicted by the isotropic elasticity. 
The core field leads to a passing angle, which depends on 
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FIG. 10. Variation of the passing stress, ay^, and of tlie pass- 
ing angle, Smax, with the dislocation dipole height, h. 

the dipole height, and which strongly deviates from 7r/4 
for small dipole heights. Such a dependence of the pass- 
ing angle with the dipole height has also been obtained by 
Henager and Hoagland for edge dislocation dipoles in fee 
metalsii^. They obtained a stronger influence of the core 
field on the dislocation interaction than in the present 
study. In their case, the core field contribution can be 
neglected only when the two dislocations are separated 
by more than 506, i.e. more than 10 nm. 



VII. CONCLUSIONS 

The approach developed in the preceding paper to 
model dislocation core fieldi^, has been applied here to 
the (111) screw dislocation in a-iron. Using ab initio cal- 
culations, we have shown that a screw dislocation creates 
a core field corresponding to a dilatation perpendicular 
to the dislocation line. The core field modeling within 
the anisotropic linear elasticity perfectly reproduces the 
atomic displacements observed in ab initio calculations. 
It also allows to derive from atomistic simulations con- 
verged values of the dislocation core energies. The de- 
veloped approach illustrates the necessity to consider the 
dislocation core field when extracting quantitative infor- 
mation from atomistic simulations of dislocations. 

Then the elastic modeling of the screw dislocation has 
been used to study the interaction energy between the 
dislocation and a carbon atom. The dislocation core field 
increases the binding of the C atom when both defects 



11 



are close enough (less than ^-^20 A). At larger separation 
distances, the C atom interacts only with the dislocation 
Volterra elastic field. 

Equilibrium properties of a screw dislocation dipole 
are also affected by the dislocation core field. This addi- 
tional elastic field increases the stability of the dipole: a 
higher stress is needed to destroy the dipole. Neverthe- 
less, when the two dislocations composing the dipole are 
sufficiently separated, one can consider that they only 
interact through their Volterra field. 

The amplitude of the dilatation corresponding to the 
dislocation core field does not depend on the dislocation 
core configuration, i.e. either easy or hard core structure. 
This indicates that the core field does not arise from 
the atomic structure of the dislocation core, but may 
be induced by anharmonicity. Our work, like previous 
similar studies^i^i^""— i^, shows that such an anharmonic 
effect can be fully considered within linear elasticity the- 
ory with the help of a localized core field. One could 
have also used non-linear elasticity theory to incorporate 
anharmonic contributions, either following the approach 
of Seeger and Haasen=l based on a Griineisen model for 
an isotropic crystal, or the iterative scheme proposed by 
Willis'^^'^^ for an anisotropic crystal. 
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FIG. 11. Core energy of the screw dislocation in the easy 
core configuration for different core radii: (a) rc = 10 A, (b) 
rc = 3 A, and (c) rc = 1 A. Both the Volterra and the core 
fields have been considered in the elastic energy. 



Appendix A: Core energies and core radius 

The core energy depends on the choice of the core ra- 
dius Tc. This core radius defines the cylindrical region 
around the dislocation line, where the strain is so high 
that elasticity theory does not apply. It therefore parti- 
tions the dislocation excess energy into two contributions, 
the core energy corresponding to the energy stored in 
this core cylinder and the elastic energy in the remaining 
space. Changing the value of modifies this partition 
between core and elastic energy without modifying the 
total excess energy [Eq. ([SJ]. In the present work, the 
choice of affects the convergence of the core energy 
with the size of the simulation box and the geometry of 
the dislocation periodic array. This arises from the dis- 
location core field. The dislocation line energy created 
by this core field depends on rc through the contribution 
1/2 MijKijkiMki 1/r-c^ [Eq. As the dipole moments 
Mij have been found to depend on the size of the sim- 
ulation cell (Fig. [3]), changing the value of rc leads to a 



shift of the core energy which also depends on this size. 
Figure [TT] shows the core energies obtained for different 
core radii. One sees that the value rc = SA leads to a 
core energy which does not depend on the geometry of 
the dislocation periodic array and which converges rea- 
sonably with the size of the simulation cell. Moreover, it 
is close to the norm of the Burgers vector (6 = 2.5 A) as 
theoretically expected^. 



Appendix B: Carbon dipolar tensor 

A solute C atom embedded in a Fe matrix is modeled 
within elasticity theory by a dipolar tensor— Pij. As 
shown in Ref. [23, the value of this tensor can be sim- 
ply deduced from the stress tensor measured in atomistic 
simulations where one solute atom is embedded in the 
solvent using periodic boundary conditions. One predicts 
that the homogeneous stress measured in these simula- 
tions varies linearly with the inverse of the volume V of 
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FIG. 12. Variation of the stress aij and of the corresponding 
dipolar tensor Pij with the inverse of the volume V o{ & unit 
cell containing one C atom in a [001] variant embedded in an 
Fe matrix. Symbols correspond to ab imtio calculations and 
solid lines to Eq. ((Bl| and (lB2l) . 



the unit cell, 



P, 



■ 

Because of the small size of the unit cell used in ab 
initio calculations, one has to take into account the 
polarizabilitj*2i"— of the C atom, i.e. the fact that the 
tensor Pij depends of the strain e'^i locally applied on the 
C atom. One can "write to first order 

"where P?- is the C dipolar tensor in an unstrained crystal 
and Plji^i its first derivative "with respect to the applied 
strain. In our atomistic simulations, this applied strain 
e'^i arises from the periodic images of the C atom. The 
strain created by a point defect varies linearly "with the 
inverse of the cube of the separation distance and we use 
for our atomistic simulations homothetic unit cells "with 
the same cubic shape. The dipolar tensor Pij should 
therefore vary linearly "with the inverse of the volume of 
the unit cell. 



P - P° 



5P^ 
V 



(Bl) 



"where SPij is a constant depending only on the shape of 
the unit cell and not on its volume. As a consequence, 
the stress measured in the atomistic simulations should 
vary like 



' V 



y2 



(B2) 



We performed ab initio calculations to obtain the val- 
ues of the dipolar tensor, characterizing one carbon atom 
in the iron matrix. We chose a cubic unit cell "which con- 
tains one C atom in an octahedral interstitial site. The 
simulation boxes contain 128, 250, or 432 Fe atoms. The 
Siesta calculation details are the same as for dislocation 



calculations, and 13 numerical pseudoatomic orbitals per 
carbon atom are used to represent the valence electrons 
as described in Ref. IH. The k-point grids used for the 
calculations are 4 x 4 x 4 for the 128 and 250 atom cells 
and 3 X 3 X 3 for the 432 atom cells. 

Because of the tetragonal symmetry of the octahedral 
interstitial site, the dipolar tensor characterizing the C 
atom is diagonal "with only t"wo independent components. 
This agrees "with the symmetry of the stress tensor given 
by ab initio calculations. The variations of this stress 
tensor "with the volume of the unit cell are in perfect 
agreement with Eq. (jB2p (Fig. [T^]). This allo"ws us to 
deduce the elastic dipolc Pij "which is characterized by 
the values Pxx — Pyy = 8.9 and Pzz — 17.5 eV for the 
[100] variant of the C atom in the dilute limit \V ^ oo 
in Eq. (IBD) ]. 

Previous ab initio calculations performed in smaller 
simulation cells have led to the C atom dipolar tensor 
deduced from Kanzaki forces^. These values are consis- 
tent "with the ones "we have deduced from the homoge- 
neous stress. 



TABLE II. Variations of the iron lattice parameter with car- 
bon atomic fraction, 5x and 5^, and formation volume of car- 
bon, 50,. 





5x 


&. 


50. (A-') 


Ah initio 


-0.086 


1.04 


10.4 


Empirical potential^^ 


-0.088 


0.56 


4.47 


Exp.ii 


-0.052 


0.76 


7.63 


Exp.i^ 


-0.0977 


0.862 


7.76 


Exp.« 


-0.09 


0.85 


7.80 



The elastic dipole, Pij, can be simply related to the pa- 
rameters 5x and 5z of the Vegard la"w^^, -which assumes 
a linear relation bet"ween the variations of the lattice 
parameter, a and c, and the atomic fraction of carbon 
atoms, xc- If all carbon atoms are located on the [001] 
variant of the octahedral site. 



a{xc) = ao (1 + 5xXc) , along the [100] or [010] axis, 
c{xc) = flo (1 + SzXc) , along the [001] axis, 

"where oq is the pure Fe lattice parameter. The param- 
eters Sx and Sz deduced from our DFT calculations are 
compared to experimental datap^"— in Table HIl The ab 
initio calculations lead to a formation volume of carbon 
larger than the experimental value, and to a larger tetrag- 
onal distortion expressed as {Sz ~ Sx). 

One can also deduce from the elastic dipole the varia- 
tion of the solution enthalpy of C in bcc Fe "with a volume 
expansion. 



Vo 



dV 



1 



- P 



V=Vo 



"where is the excess enthalpy of a Fe crystal of equi- 
librium volume Vq containing one C atom. The value 
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corresponding to the elastic dipoles calculated above, 
— ll.SeV, is in good agreement with the value obtained 
by Hristova et al^, -12.3 eV, using an ab initio ap- 
proach based on GGA-DFT with a plane waves basis set 



and the Blochl projector-augmented wave method (PAW) 
as implemented in the Vienna Ab Initio Simulation Pack- 
age (VASP). This validates our Siesta approach for char- 
acterizing C atom embedded in an iron bcc matrix. 
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